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ABSTRACT 


Detailed measurements of the axial velocity profile and 
electromagnetic structure of a high power, quasi-steady MPD 
discharge are used to formulate a gasdynamic model of the 
acceleration process. Conceptually dividing the accelerated 
plasma into an inner flow and an outer flow, it is found that 
more than two-thirds of the total power in the plasma is de- 
posited in the inner flow, accelerating it to an exhaust veloc- 
ity of 12.5 km/sec. The outer flow, which is accelerated to a 
velocity of only 6.2 km/sec, appears to provide a current con- 
duction path between the inner flow and the anode. 

Related cathode studies have shown that the critical cur- 
rent for the onset of terminal voltage fluctuations, which was 
recently shown to be a function of the cathode area, appears 
to reach an asymptote for cathodes of very large surface area. 
Detailed floating potential measurements show that the fluctua- 
tions are confined to the vicinity of the catnode and hence re- 
flect a cathode emission process rather than a fundamental limit 
on MPD performance. 

A 45-cm wide parallel plate configuration produces a two- 
dimensional discharge which replicates the essential features 
of the conventional three-dimensional MPD discharge. Magnetic 
probe studies show that for a total current of 65 kA, the cur- 
rent distribution is uniform i-o better than 10% across the entire 
electrode width with no evidence of instabilities. Using a re- 
cently constructed resonant optical cavity which will straddle 
this configuration and a newly designed and calibrated mass feed 
system capable of injecting a uniformly distributed mass flow of 
more than 10 g/sec with a risetime of only 2 msec, this device 
will be capable of determining directly whether lasing can be 
obtained in the MPD flow. 
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Axial floating potential and enclosed current distribu- 
tions have been measured for several hollow cathode configura- 
tions using an extended current pulse of 500 psec duration. 

For a cathode cavity which diverges with a 10° half-angle, 
the characteristic features of the distributions are the same 

as for previous cylindrical hollow cathodes: a surface current 

2 

density greater than 1000 A/cm at the downstream end of the 
cavity accompanied by a weak axial electric field of less than 
10 V/cm. For a configuration where the cavity diverges at a 
45° half-angle, the equipotential plateau recedes further in- 
side the cavity and drops to only +4 volts with respect to the 
cathode, as compared to approximately +20 volts for the cylin- 
drical cathode. In addition, preliminary data indicate that the 
current density peaks on the 45° face rather than at the down- 
stream end of the cavity. Based on the results of these and 
previous tests, a cylindrical hollow cathode with refractory in- 
sulation has been selected for a more detailed series of diagnos- 
tic measurements. 
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I. INTRODUCTION 

During this past spring, the Electric Propulsion Lab 
was moved from the Forrestal Campus, approximately four miles 
outside of Princeton, to the Engineering Quadrangle on the 
main campus. The move was undertaken to facilitate collabora- 
tion with several related research programs already housed in 
the Quadrangle, and to provide more accessible facilities for 
student participation in independent research programs during 
their undergraduate years. _ The relocation was totally financed 
by the University, which continues to regard the Princeton pro- 
gram in electric propulsion as one of the principal vehicles 
for undergraduate and graduate research education within the 
Engineering School, and one of the selected areas of engineer- 
ing research in which it can best serve the national aerospace 
needs. The new laboratory was officially opened on Tuesday, 

May 21 in a brief ceremony attended by the President of the 
University, uilliam Bowen, Mr. James Lazar from NASA Headquarters, 
and approximately 50 members of the Engineering faculty, staff 
and students. Figure 1 shows Mr. Lazar triggering the initial 
plasma discharge in the new laboratory. 

During Mr. Lazar's visit, the various phases of the labora- 
tory research program were reviewed in detail. After discussing 
the on-going studies of power deposition into the MPD anode and 
the insulator ablation and cathode emission problems, the sug- 
gestion was made that the MPD discharge work be broadened to 
serve as a basis for further arc physics studies. Accordingly, 
the MPD effort now encompasses the Ph.D. programs of three grad- 
uate students and the independent research topic of one under- 
graduate student. Of these, the work on the acceleration mechan- 
isms in an MPD discharge free of ablated species, and two aspects 
of the study of the cathode .emission processes are reviewed in 
this semi-annual report. 
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FIRST FIRING OF DISCHARGE APPARATUS BY 
J. LAZAR, NASA HEADQUARTERS 


FIGURE I 

AP 25 P-539 
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The plasmadynamic laser program is an example of the de- 
velopment of a new avenue of research from a more general otudy 
of the physics of high current MPD discharges. Earlier measure- 
ments of the radiation environment in the discharge revealed 
several levels of ionization in the expanding exhaust flow. 

Use of operating conditions where this radiation is optimized 
gives rise to the possibility of high power lasing in the 
visible and ultraviolet. Of the two Ph.D. programs in this 
topical area, one is centered on the examination of the optical 
depth of selected argon ion and neutral lines in the exhaust 
of the conventional three-dimensional discharge geometry. The 
second program examines the gain of selected lines directly by 
using an optical cavity in conjunction with a unique two-dimen- 
sional plasma accelerator geometry. This latter program has 
progressed rapidly and will be described in detail. In addition 
and at the sponsor's request, technical contact in this area has 
been initiated with F. C. Schwenk and K. Thom of NASA Headquarters. 

The program to determine the physics of hollow cathode opera- 
tion has completed its first phase. This work, which is reviewed 
herein, has exariined various hollow cathode geometries and the 
extent to which they influence the current distribution within 
the cavity and the equipotential profile along the cathode center- 
line. Based on these results, a cathode has been selected with 
which to begin a series of more detailed measurements of the plasma 
properties within the cavity. 
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I I . LABORATORY RELOCATION 

Following the recent move, the Electric Propulsion Labora- 
tory is now located on the main campus in the Engineering Quad- 
rangle, shown in Fig. 2. This building houses the classrooms, 
experimental and support facilities for the Electrical Engineer- 
ing, Chemical Engineering, Civil Engineering, and Aerospace and 
Mechanical Sciences Departments, as well as several interdis- 
ciplinary programs in energy conversion, bioengineering, environ- 
mental studies and transportation. 

The move proved to be an advantageous time to redesign the 
laboratory layout for optimum service of the present and future 
needs and for incorporating several modifications and improve- 
ments to the existing discharge equipment. Figures 3, 4 and 5 
show three views of the new laboratory arrangement. As can be 
seen, large areas are now available around each machine for 
improved accessibility and specialized diagnostic equipment. 

In addition, sufficient flexibility has been built in to allow 
repositioning each machine within its area for special tests 
or discharge configurations. The photos also show an optics 
room connected to the main lab and a conference area, two 
features that were not present in the previous laboratory 
arrangement due to space limitations. Not shown but immediate- 
ly adjoining the optics room is a photographic dark room. 

All of the capacitor power supplies are located on the 
floor below the laboratory with the connection to the dis- 
charge devices provided by bundles of coaxial cable (RG 8/U) . 
Figure 6 shows the two principal capacitor banks. The bank 
on the left provides currents from 1 to 20 kA for hollow cath- 
ode experiments, and the bank on the right supplies up to 
120 kA for MPD discharges associated with either propulsion, arc 
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VIEW OF HOLLOW CATHODE FACILITY 




FIGURE 6 
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physics or plasmadynamic laser studies. Both capacitor banks 

and all. experimental facilities in the main lab are now fully I 

operational following this relocation. j 
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III. QUASI-STEADY MPD DISCHARGE 

A. Plasma Acceleration Processes (Boyle) 

In order to determine the maximum performance capability 
of the MPD accelerator, it is necessary to understand com- 
pletely the plasma acceleration mechanisms within the discharge. 
However, previous semi-annual reports have shown that the ef- 
fects of insulator ablation and cathode surface area can great- 

i 5 ]^ ]_55 

ly alter arc behavior and thus obscure this evaluation. * ' 

With these effects reduced to a minimum, it is now possible to 
proceed systematically toward acquiring this understanding. 

In principle, the total velocity vector field, and thus 
the acceleration pattern, of the MPD device may be determined 
by the simultaneous solution of the mass, momentum and energy 
conservation equations for given initial and boundary conditions. 
The electromagnetic variables responsible for the Lorentz body 
force in the momentum equation and the Joule heating term in 
the energy equation must be consistent with Maxwell's equations. 
In addition, a generalized Ohm's lav; which describes the coupling 
between the electromagnetic and gasdynamic variables must be 
prescribed. An appropriate equation of state and expressions 
for the plasma transport properties including a plasma con- 
ductivity law are required to complete the set. In practice 
the coupled nonlinear nature of the governing partial differ- 
ential equations presents a proble/r. of extreme mathematical 
complexity. 

By contrast, the quasi-steady plasma acceleration problem 
becomes more tractable when discussed in a semi-empirical 
manner. In this section the measured exhaust velocity profile 
and the electromagnetic discharge structure are presented for 
an accelerator operating at an arc current of 15.3 kA and an 
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argon mass flow of 6 g/sec. These data are combined in a 
quasi-phenomenological model of the acceleration process and 
finally discussed in relation to the fundamental capability 
of the quasi-steady accelerator. 

Experimental Conditions 

The plasma acceleration process may be studied for any 

desired combination of arcjet parameters so long as the 

accelerator operation is uninfluenced by ablation or cathode 

surface area effects. There is particular interest, however, 

in examining the exhaust velocity profile of the quasi-steady 

accelerator for operating conditions such that the parameter 
2 

J /m, the arc current squared divided by the argon propellant 

mass flow rate, equals approximately 40 kA -sec/g. Claims 

have been made previously that the acceleration mechanisms are 

A— 1 

fundamentally limited at these conditions. Furthermore, re- 

sults in the last semi-annual report indicate that arcjet 

operation free from (boron nitride insulator) ablation is re- 

2 2 
stricted to J /m less than or equal to 40 kA -sec/g. 

This condition is realized experimentally with an arc 
current or 15.3 kA and a total argon mass flow rate of 6 g/sec. 
The mass injection configuration equally divides the total 
propellant flow between inner and outer injection orifices. 

Three grams per second are injected directly about the cathode 
through an inner cathode base annulus. The remaining 3 g/sec 
enter the discharge chamber through twelve 0.32 cm holes sym- 
metrically distributed about the 3.81 cm radius. 

To prevent voltage fluctuations, a cathode area greater 
2 

than 20 cm is required for this combination of arc current and 
155 

mass flow. A thoriated tungsten cathode 1.91 cm in diameter, 

7.62 cm in length, with a 2.54 cm conical tip provides sufficient 

2 

electrode area («41 cm ) to insure stable arcjet operation. 

The terminal arc voltage for this configuration equals 115 
volts. Typical current and voltage signatures at these con- 
ditions are shown in Fig. 7. 



ARC CURRENT AND VOLTAGE SIGNATURES 


FIGURE 7 

AP25-P-526 
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Axial Velocity Profile 

The axial velocity profile measured using the time-of- 
flight probe technique on the accelerator's centerline is 
presented in Fig. 8. The axial distance z is measured 
downstream from the anode plane. Velocity measurements 
within the discharge chamber were precluded due to the lack of 
discernable ion density fluctuations required by the time-of- 
f light technique. The plasma leaves the discharge chairiber at 
a velocity less than 9.0 km/sec. Final exhaust velocities of 
12.5 km/sec are observed several anode orifice diameters down- 
stream. This final exhaust velocity exceeds by nearly 50 % the 
Alfven critical speed of 8.7 km/sec for argon, previously sug- 
gested as the upper limit to which propellant may be accelera- 
ted by a quasi-steady MPD arc operating at these conditions. 

Electroma gnet ic Stru cture of the MPD Discharge 

The exhaust velocity profile results primarily from the 
combined action of electromagnetic body forces upon the accel- 
erated propellant mass. The magnitude and direction of the 
electromagnetic body forces may be ascertained from measured 
distributions of current density and magnetic field. In this 
way the local acceleration processes may be related to the 
local j x B distribution. Alternatively, the plasma accelera- 
tion process may be viewed as the conversion of some portion 
of the electrical input power (j’E) into directed motion of the 
propellant. Evaluation of the local j • E power deposition within 
the arc discharge requires knowledge of the electric fields as 
well as the current density distribution. From these distribu- 
tions of magnetic field, current density and electric field, a 
qualitative model of the quasi-steady acceleration mechanisms 
will emerge. 



VELOCITY, km/sec 
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CENTERLINE AXIAL VELOCITY PROFILE 


FIGURE 8 

AP25-5050 
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Because the arc discharge is cylindrically symmetric, the 

* 

induced self magnetic field (B) is azimuthal. The B field 
distribution is determined from axial and radial induction 
coil surveys. Azimuthal probe surveys for selected radial 
locations within the discharge confirm the presumed cylin- 
drical symmetry of the magnetic field. The magnetic field 
distribution stabilizes in two characteristic times. Approxi- 
mately 200 psec after discharge initiation, the magnetic field 
within the discharge chamber assumes a quasi-steady distribution. 
At a somewhat later time, approximately 450 psec into the 
1000 psec current pulse, a quasi-steady distribution throughout 
the entire flow field is achieved. This quasi— steady distribu- 
tion, presented in Fig. 9, is maintained for the duration of 
the current pulse. In this figure, lines of constant magnetic 
induction are superimposed upon a schematic of the accelerator 
chamber . 

The quasi— steady current distribution is derived directly 
from the magnetic field data employing Ampere's law and the 
azimuthal symmetry of the discharge. The resulting distribution 
of enclosed current for the total discharge current of 15.3 kA 
is illustrated in Fig. 10. Each current contour represents a 
current streamline and is labelled with the amount of current 
in kiloamperes which flows downstream of it. The magnitude of 
the local current density vector may be estimated by dividing 
the current between two adjacent current contours by the normal 
area through which it flows, 

A significant fraction of the total arc current flows in 
the region exterior to the discharge chamber. The current 
streamlines within the accelerator chamber emerge from the cath- 
ode with a density distribuf-ion peaked at both ends of the cath- 
ode as shown in Fig. 11. The current contours then bend axially 
downstream be 'ore turning radially toward the anode. Anode cur- 
rent attachment concentrates about the lip region of the anode. 
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Outside the discharge chamber nearly fifty percent of the 
total arc current billows far downstream in the exhaust plume 
before returning and attaching on the outer surfaces of the 
anode barrel. The current density vector within several centi- 
meters of the accelerator' s centerline in this exhaust region 
is primarily axial. This current density pattern is charac- 
teristic of the moderate magnetic Reynolds number, 6(1-10), 
associated with the discharge flow. 

The local electric fields, which integrate to the terminal 
arc voltage, arise as a result of current conduction across the 
flowing, finite conductivity discharge plasma. Current conduc- 
tion processes across the electrode-plasma interface also create 
potential gradients commonly referred to as electrode falls. 

The local electric field equals the gradient in plasma 
potential which, in turn, is related to the experimentally 
measurable floating potential. For small distances over which 
the gradients in electron temperature, flow vector, and col- 
lision frequencies are negligible, the electric field may be 
evaluated from the difference in floating potential. The quasi- 
steady distribution of floating potential contours is shown in 
Fig. 12. The cathode is at 115 volts negative with respect to 
the anode ground. The labelled equipotential contours are also 
negative with respect to the ground, the negative signs being 
omitted from Fig. 12 for convenience. The 20 and 90 volt con- 
tours envelope the anode and cathode regions respectively. The 
equipotential contours in the outer radial regions of the down- 
stream exhaust are approximately parallel to the current stream- 
lines. 

Evaluation of the total arc jet power distribution through- 
out the MPD accelerator requires local evaluation of j"*E, the 
dot product between the current density and electric field vec- 
tors. For this purpose, the accelerator may be sub-divided in- 
to several separate regions, A through F, as is done in Fig. 13. 




0 5 10 15 

Z) CR1 

QUASI - STEADY FLOATING POTENTIAL CONTOURS 


05^3 

gm 

ro 


r,cm 




mi 





1 

ro 


N 

POWER 

DEPOSTION ' 

REGION 

POWER, watts 

A 

6.IXI0 3 

B 

2.5XIO® 

C 

2.0X10° 

D 

— 

E 

4.4X10° 

F 

4.4X10° 




23 


These regions, excluding the electrode fall regions E and F, 
are further divided into incremental volume elements in which 
average values for the magnitude and direction of the j and E 
vectors are assigned. The product of local power density and 
incremental volume element when summed over each region repre- 
sen. s the total electrical power deposited in that region. 

This regional power deposition is summarized in Table I . At 
15.3 kA and 6 g/sec the total arc power is 1.76 megawatts. 

Table I 


Region 

A 

B 

C 

D 

E 

F 

Total Arc Power 


(5-E) (vol) , watts 

6.1 x 10 5 
2.5 x 10 5 
2.0 x 10 4 

4.4 x 10 5 
4.4 x 10 5 
1.76 x 10° 


Region description 

inner chamber flow 
outer chamber flow 
inner exhaust flow 
outer exhaust flow 
anode fall 
cathode fall 


To first order, no power addition appears in region D of 
the discharge because the electric field is perpendicular to 
the current density vector, thereby making j-E zero throughout 
the volume. The power lost to the anode in region E is estimated 
from recent work relating anode power fraction to total arc power 
in the quasi-steady MPD discharge. ^ ^ For J^/m = 40 kA^-see/g, 
these results extrapolate to an anode power fraction of 25%; i.e. 
25% of the total arc power or 4.4 x 10 watts is accounted for 
by loss mechanisms in the vicinity of the anode. With the anode 
fall and bulk plasma flow regions accounted for, the remaining 
input power must appear in the cathode fall region F to provide 
an overall power balance of the arc. 
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The distribution of local power deposition for the MPD 
accelerator operating at 15.3 kA and 6 g/ sec may now be 
summarized: Approximately fifty percent of the total arc power 

is invested in the plasma flow which occupies regions A, B, 
and C in Fig. 13. The major power addition to the plasma flow 
occurs within the discharge chamber, with negligible power 
addition ensuing in the downstream exhaust. The remaining 
fifty percent of the total power appears in the electrode fall 
regions of the discharge. 

Flow F ield Characteristics of the MPD Discharge 

The flow pattern in the exhaust of the accelerator is 
determined from velocity vector measurements using the double 

g c 

probe. The mass streamlines constructed from these local 

flow angles parallel the current streamlines in the inner ex- 
haust region C depicted in Fig. 13. Precise flow vector measure- 
ments in the outer exhaust region D are difficult to obtain due 
to a decrease in the probe's sensitivity to the fl-.w orientation. 
Despite this decreased accuracy, flow vectors not parallel to the 
current streamlines are identified. The mass streamlines appear 
to diverge more rapidly than the current streamlines in the outer 
exhaust region. 

The flow in region C of the exhaust is supersonic. Bow 
shocks are photographed off the hemispherical tips of rods intro- 
duced into this inner flow several centimeters downstream of the 
anode orifice. Oblique shocks attached to wedges of different 
turning angles are also observed. As these wedges are moved 
axially inward towards the anode plane, the oblique shock angles 
increase indicative of lower incident flow Mach numbers. Finally 
no shocks are seen in the vicinity of the anode orifice indicating 
near sonic conditions there. 
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Insufficient knowledge of the shocked plasma's thermo- 
dynamic state precludes exact Mach number calculations from 
the shock angle data. However, the oblique shock relations 

for a perfect gas with an appropriate ratio of specific heats 

A— 2 

provide a first order estimate. Using this method, the 

flow Mach number seven centimeters downstream of the anode 
orifice is estimated between two and three. 

Photographs through narrow band spectral filters display 

the luminosity patterns of the discharge. Spectrograms of the 

discharge reveal singly ionized argon (A-Il) as the predominant 

species. Argon II line radiation in the downstream exhaust 

extends to radii beyond the anode barrel of the accelerator 

although the spectral intensity monotonically decreases with 

radius. Referring again to Fig. 13, the predominant argon 

o 

II radiation photographed through a 4880 A filter is confined 
to the inner flow regions A and C. The outer edge of the A-II 
luminosity in region C coincides with a mass streamline pre- 
viously defined by the velocity vector measurements. The cor- 

o 

respondence among the 4880 A A-II luminosity boundary, the mass 
streamline and the demarcation of the A-C and B-D flow regions 
is illustrated in Fig. 14. 

The remaining regions of the discharge, B and D, are 
essentially dark, except for a small toroidal volume located 
in the anode plane between the anode lip and region A. 

Quasi-Ste a dy MPD Acceleration Processes 

Having discussed the pertinent electromagnetic and flow 
field characteristics of the MPD discharge, a model of the 
acceleration mechanisms responsible for the measured plasma 
velocity profile may be constructed. A complete description 
of the flow processes would require additional profiles of 
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number density, ion and electron temperatures, and pressure 
throughout the discharge. For this arc jet configuration, 
these profiles have not yet been determined. However, 
reasonable estimates of these variables based upon previous 
measurements for other configurations can be made. This 
allows a phenomenological model of the plasma acceleration 
process to be formulated. 

Two— Flow P la sma Acceleration Model 

The total injected argon propellant is divided into two 
independent mass flows associated with chamber injection at 
large and small radii. The local power deposition and lum- 
inosity profiles in regions A and B of the discharge chamber 
suggest this flow division is essentially maintained as the 
two flows are accelerated axially downstream out of the cham- 
ber and into the exhaust regions. The plasma acceleration 
process will thus be modelled in terms of two separate inner 
and outer propellant flows, schematically illustrated in Fig. 

15. 

Inner and outer flow separation is consistent with the ob- 
served electromagnetic discharge structure within the accelerator 
chamber. In particular, there are substantial differences 
between the axial power deposition profiles for the two flow 
regions. Figure 16 compares the incremental power addition 
to each flow as a function of axial distance downstream from 
their respective injection orifices. Approximately 180 kW are 
deposited immediately in the inner flow as it emerges from the 
cathode base injection annulus. This power level is commen- 
surate with complete single ionization of the injected 3 g/sec 
argon. Once ionized, the prevailing j x B body forces tend 
to confine radially and accelerate axially the inner flow. 

Only after leaving the discharge chamber does the inner flow 
expand radially. 







AP t , watts 



AXIAL DISTANCE, cm 

INCREMENTAL INPUT POWER PROFILES 

FIGURE 16 

AP25»5054 



30 


In contrast, negligible input power appears in the outer 
flow upstream of the anode lip region. A total of only 27 kW 
is deposited in this region. Thus, it may be inferred that 
only a small fraction of the 3 g/sec outer mass flow rate is 
ionized upstream of the anode. This is also consistent with 
the absence of ionized argon line radiation in this region of 
the discharge. In addition the local current density and mag- 
netic field are substantially lower in this region. As a con- 
sequence of the weak local Lorentz forces, the initial axial 
flow direction of the outer propellant stream will be only 
slightly' perturbed. To this approximation, the integrity of 
the two concentric inner and outer flows is preserved as they 
are accelerated downstream. 

Inner Flow Region . Two— thirds of the electrical input power 
invested in all the plasma flow regions appears as total enthalpy 
of the inner propellant flow. In addition, nearly 100 percent of 
this inner flow power addition is accomplished within the dis- 
charge chamber as shown in Fig. 17. Negligible power addition 
occurs in the exhaust. The discharge chamber and exhaust plume 
regions of the inner flow are distinguishable by this localiza- 
tion of power input. As such, the flow processes within each 
will be treated separately. 

Inner Flow, Discharge Chamber . The inner propellant mass 
flow is ionized, heated, and accelerated through the sonic point 
within the discharge chamber. The net result of these processes 
may be examined in terms of a power balance across the region. 

The total input power, P , invested in the inner flow with- 

m the discharge chamber is 6.1 x 10 watts, exclusive of the 

cathode fall power. Of this total, a certain amount, P , is 

R 

lost through radiative processes. The remaining power appears 
directly in the flow. 
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P FLOW ~ Pr 

This flow power is partitioned between enthalpic and 
kinetic modes. 

P*o„ = ™(h+ U Z /2) 

where ft is the accelerated mass flow rate, h the flow 
enthalpy, and u the mean flow velocity. 

A quasi— equilibrium expression for the flow enthalpy 
is assumed, in the form, 

h = fc p (r,p)Jr +JLRfjT 

o % 

where T g and T^ are the electron and ion temperatures, 

Cp is the heat capacity at constant pressure and R equals 
the argon gas constant. The first integral in Eqn. III-3 
represents the contribution to the flow enthalpy from both 

excitation and ionization as well as the energy content of 

A— 3 

the ion and electron random motion at temperature T . 

Generally the excitation and ionization processes in devices of 
this type limit the electron temperature to between one and 
two volts. After initially baring ionized, the plasma exper- 
iences a significant axial acceleration. The resulting high 
velocities are characterized by flow times short compared to 
characteri st ic electron— ion energy equipartition times^ - ^ 
thereby establishing a non- -equi librium situation. Because 
electron— ion collisions are unable to establish an squi 1 : briv >.i 
temperature distribution in time, additional enthalpy appears 
in ion random motion alone, thereby raising the ion temperature 
above the electron temperature. This non-equilibrium effect 
is accounted for by the second integral in Eqn. III--3. 


(III-l) 


(III-2) 
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Equations III-l, 2 and 3 may "be rearranged to express 
the flow velocity in terms of input and radiated power, 
accelerated mass flow rate, pressure, and electron and ion 
temperatures such that: 


u= \z( Pt-p« + sR(r e -T L ) ~fc p lr,p)dr^ 

l \ rij z J 0 


j] 




(III-4) 


Order of magnitude estimates of P , p, T and T. at a 
given axial location are required to calculate the local flow 
velocity from the measured power input up to that point. The 
measurement techniques for any one of these variables require 
a major experimental effort in terms of the sophistication 
and time demanded by the particular measurement. In lieu of 
their actual measurement, these estimates' must rely upon pre- 
vious data obtained at similar operating conditions, albeit 
for different arc jet configurations. 


Plasma radiation phenomena in an environment of this 
type are extremely complex. An exact formulation of the in- 
put power radiated away by bremmstrahlung and line radiation 
is beyond the scope of this work. The power density radiated 
in free— free transitions for a Maxwellian electron distribution 
is given by Spitzer as: 


P H (w/cn 3 ) = /.4Z(to' 34 )z z n e ni%' A (HI-5) 

where Z is the ion charge number and n . are the electron 

A-4 e ' x 

and ion number densities. For a singly ionized argon plasma 

15 — 3 

with an electron number density equal to 10 cm and T g be- 
tween one and two volts: 

P H = 0(1 HI /cm 3 ) 

2 3 

The power radiated away by bremmstrahlung from the 10 cm 
inner flow volume may thus be approximated as: 


■••'‘ISO* 
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p" = P„ (</ol) = 0 (to 2 w) 

a negligible amount relative to the total input power. How- 
ever, power losses from bound-bound transitions may be sig- 

2 

nificant. A line radiation flux of 2000 watts/cm has pre- 

125 

viously been found typical of this plasma environment. 

The luminous inner flow region within the discharge chamber 

2 2 

radiates from a surface area of 0(10 cm ). From this sur- 
face , 

P R b * 0 (zx / 0 s jv) 

will be radiated away. The assumption of no recombination 
within this discharge region implies the power loss from free- 
bound transitions is negligible. The total input power lost 
to radiation is therefore 

P R = P R ff -h P R b = 0(z*io s vs) 

Ele<- tron temperatures consistent with previous probe 

measurements and observed levels of ionization, pressures 

inferred from static pressure measurements in the discharge 

chamber, and argon ion temperatures extrapolated from Doppler 
. . . ...... 125,134,141 

width data are listed below. 

T e = 0(19,000 °K) 
p = 0 (10 _2 -10 _1 atm) 

T i = 0(60,000 °K) 

5 

For P T = 6.1 x 10 watts just outside the discharge chamber 
and an accelerated inner mass flow rate of 3 g/sec, the cor- 
responding exit velocity is calculated from Eqn. Ill— 4 to be 
8.0 ~u — 10.4 km/sec. The experimentally measured plasma 
velocity 2 centimeters downstream is 9.6 - 0.4 km/sec. The 
inner propellant flow thus is' slightly supersonic as it de- 
parts the discharge chamber with an elevated ion temperature 






jii'. 
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and an electron temperature consistent with the conductivity 
requirements of the discharge. 

Inner Flow, Exhaust Plume . The flow in the inner exhaust 
plume of the discharge is accelerated supersonically to a final 
velocity of 12.5 km/sec. This acceleration region extends 
several anode orifice diameters downstream. Within this region 
the electromagnetic body forces are normal to the flow stream- 
lines. Thus, electromagnetic body forces tend to conf : ne 
rather than accelerate the flow in the streamwise direction. 

The axial acceleration, however, may be accounted for by the 
radial divergence of the not, electromagnetically confined 
inner plasma flow. 

The mathematical description of the expansion process is 
significantly simplified for one-dimensional flows. Radial 
profiles of velocity and number density at several axial loca- 
tions justify the 1-D assumption as a first approximation of 

the flow field. The one-dimensional mass, momentum, and energy 

A— 5 

relations for a given streamtube may be written as: 


J( f> u A) - 0 
j>u da + dp = 0 
dh + u da = dc 


(III— 6) 
(III— 7) 

(III-8) 


where A is the normal streamtube area, dq the energy addition 
per unit mass from Joule heating, and , u, p, and h the common 
nomenclature for density, velocity, pressure, and enthalpy. The 
flow enthalpy may be expressed as 


h = | P/f + h i 


(III-9) 


h^. is the internal energy content due to excitation and 


where 

j. 

ionization. The energy addition, dq, to the expanding flow was 
found to be insignificant compared to the total energy content. 
Thus, dq may be approximated as zero. 
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Significant argon ion recombination (AIII -'•All or 

AII-^AI) in this exhaust region of the discharge has not 

A— 6 

been experimentally observed. Recombination times for 

these processes are characteristically long, compared to 
local flow times. This suggests the plasma flow may be 
treated as frozen, i.e., 

dh L - 0 

Hence, the one-dimensional energy equation becomes 

d($ U du- ** 0 (in-10) 

Equation III— 10 integrates to 

5 p/p + u z /z =* K (iii-ii) 

2 f 


where the integration constant K represents the total 
specific flow enthalpy exclusive of excitation and ionization 

energies. The constant K may be evaluated from either the 

* 

sonic velocity u or the maximum velocity u according to 

max 

Eqn. Ill— 12 . 

Z * 


K = 5+$ u * 
2* 


= U 


rn&x 


(III-12) 


where the sonic condition 


Yp*/f* = u 


*2 


has been used. 


Equations 111-6,7, 10 and 12 may be solved for the stream- 

★ 

line area as a function of u, K, and the minimum area A . It 
is easily shown that 


A_ 

A* 


K - 


5 /2 


K - uVz 

or, in terms of the streamtube radius r, 

r * r * [*- u*Vz 1 /4 

K - « 2 /k J 



(III— 13) 



(III— 14) 
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Mass streamline coordinates r(z) for given. r* values may be 

calculated when the measured axial velocity profile u(z) in 

Fig. 8 is substituted into Eqn. 111-14. K and u* are evaluated 

from u =12.5 km/sec. The flow streamlines associated with 
max 

r* = 1.0, 1.5, and 2.0 centimeters are shown in Fig. 13, Cir- 
cular and triangular data symbols represent two experimentally 
measured streamlines. The circular data symbols also represent 

the outermost boundary of the All luminosity photographed through 
o 

the 4880 A filter. Streamline agreement between the one- 
dimensional expansion model and the experimental data is excel- 
lent. Thus, the downstream acceleration of the inner propellant 
flow may be accounted for by the radial expansion of the high 
enthalpy plasma discharged from the accelerator chamber. 

Questions may arise as to the sensitivity of these results 
to the shape of the measured velocity profile and the assumption 
of negligible recombination in the expanding flow. The above 
calculations were repeated for an arbitrarily chosen linear 
velocity profile substituted for the measured profile of Fig. 8. 
The corresponding streamline shapes calculated from Eqn III-14 
for the linear profile differed noticeably from the experimental- 
ly measured flow directions. Recombination enthalpy was arbi- 
trarily introduced onto the flow by including dh^ in Eqn. III-10. 
Even ^or relatively small degrees of assumed recombination, the 
calculated streamlines bore no resemblance to the measured flow 
patterns . 

O uter Flow Region . Plasma acceleration appears qualitative- 
ly less significant in the outer flow regions of the discharge. 
Reduced levels of power addition, local j x B, and luminosity 
suggest lower velocities, densities, and ion temperatures in 
the outer flow. Estimates of these temperatures and velocities 
may be made in the same spirit as with the inner flow. 

Outer Flow, Discharge Chamlier . The interesting feature of 
the outer flow within the discharge chamber is the localization 


I 
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of power addition off the tip of the anode previously il- 
lustrated in Fig. 16. Because of this, significant ioniza- 
tion of the injected propellant is delayed until the flow 
enters the vicinity of the anode orifice. The local power 
addition there is commensurate with full ionization of the 
incident propellant flux. It should be recalled that outer 
flow mass injection was originally required to prevent severe 
terminal voltage fluctuations and minimize the total arc 
voltage. Similar behavior has been reported for different 
arcjet configurations operating at conditions where the anode 
region appears locally "starved" of propellant. One may 
speculate that the outer flow functions primarily to provide 
a source of conduction electrons in the vicinity of the plasma- 
anode interface. 


The outer flow plasma velocity in the anode plane may be 
estimated using the power balance Eqn. III-4. However, ion 
temperatures in this region of the discharge have not been 
previously measured. To circumvent this problem, advantage 
may be made of the near sonic conditions existing in the 
vicinity of the anode orifice. The ion temperature may be 
expressed in terms of the acoustic speed c s . 


c/ ‘R (*e T e + *,• rj 


where Y • are the respective electron and ion specific 
heat ratios. Since the acoustic speed equals the flow 
velocity at the sonic point, the ion temperature may be 
expressed as. 


T i - 


U 

r ; R 


- — r e 


(III-15) 


Using Eqn. Ill— 15, Eqn. III-4 may be rearranged to express 
the flow velocity u in terms of P„ T< _ T ,, m, and T evaluate 

r JjOW 0 

at the sonic point. Electron temperatures on the same order 


1 
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as the inner flow are assumed, i.e. 0(19,000 °K) . The elec- 
trical power input to the outer flow within the discharge 
chamber is 2.5 x 10 5 watts. This input power level according- 
ly accelerates the 3 g/sec outer flow rate to velocities less 
than 4500 m/sec and heats the ions to temperatures below 
40,000 °K. Both these values are approximately a factor of 
two less than the inner flow velocity and ion temperature at 
a similar axial station. 

Oute r Fl ow, Exhaust Plume . The mass streamlines measured 
in the outer exhaust plume flow of the accelerator do not paral- 
lel the local current streamlines. The more rapid radial di- 
vergence of the outer exhaust flow indicates only partial mag- 
netic confinement. The expansion process in this region of the 
discharge should therefore more closely resemble that of a free 
expansion from a sonic orifice. The details of the outer ex- 
haust flow field were not considered. 

Specific Impulse of the MPD Discharge Flow 

The specific impulse of the quasi-steady plasma accelerator, 
defined by Eqn. III-16, is an important parameter used to evalu- 
ate thruster performance. 

% sp ~ T / (III— 16) 

where T is the thrust produced by the acceleration of the 

mass flow rate m, and g is the standard acceleration of 

2 

gravity. The AVCO critical J /m model implies an upper limit 

of 890 seconds for the quasi-steady acceleration of argon 

A— 1 

propellant. The foregoing results of the two-flow plasma 

acceleration model m*y be examined in light of this claim. In 
order to do this the thrust producing regions of the discharge 
must first be identified. 

The thrust produced by the quasi-steady MPD arcjet is both 
electromagnetic and electrothermal in origin. Electromagnetic 


! 


thrust arises from vhe coupling among the vector components of 
the current density and azimuthal self-magnetic field of the 
discharge. Both axial and radial Lorentz body forces result 
from radial and axial current flow espectively . The 3 r ^Q 
body forces ("electromagnetic blowing") directly accelerate 
the propellant downstream imparting thrust to the arcjet 
through the interaction between the current flow in the dis- 
charge plasma and the accelerator circuitry. The magnitude 
of the "blowing" thrust component may be calculated by inte- 
grating the Maxwell stress tensor over any convenient surface 

6 1 

which encloses the discharge current pattern. The radial 

or "pumping" j B„ forces confine the discharge plasma and 
support radial pressure gradients similar to those found in 
classical pinch discharges. The thrust contribution from 
these "electromagnetic pumping" forces may be evaluated by 
integrating the pressure distributions caused by this magnetic 
confinement over those surfaces normal to the accelerator’s 
symmetry axis. Thrust results from the overpressure generated 
on these surfaces. 

Electrothermal or aerodynamic thrust generation is a 
consequence of the discharge plasma being resistively heated. 

As a result, the unbalanced pressure on the discharge chamber 
walls and any subsequent nozzle-like flow expansion in the 
exhaust, vill impart additional thrust to the device. The 
proper evaluation of the electrothermal thrust component is not 
straightforward. Conventionally, aerodynamic thrust is ex- 
pressed in terms of a nozzle discharge coefficient. However, 
the appropriate discharge coefficients and chamber pressures 
for the quasi-steady MPD arc are generally unknown. Usually 
therefore, one roust rely upon specific experimental results 
to estimate the magnitude of the electrothermal thrust com- 
ponent . 
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The thrust attributable to the inner and outer flows may 

be estimated from the experimental data taken at J = 15.3 kA 

1 58 

and m = 6 g/sec. The results are summarized below. 

Inner flow: 


T 1 

= 

T 1 

em 

+ 

T 1 

et 


• ~* 

28.0 

+ 

9.5 



37.5 

N 



Outer flow: 


T° 


T° 

em 

+ 

T°, 

et 


= 

16.0 

+ 

2.7 


— 

18.7 

N 



where the superscripts I and 0 denote the inner and outer flows 
and the subscripts em and et imply electromagnetic and electro- 
thermal thrust components. The total thrust generated by the 
MPD discharge at these conditions is the sum of T -1 and T u and 
equals , 

T = T 1 + T° 

= 56.2 N 


It is interesting to note that approximately 22% of this total 

calculated thrust is electrothermal in origin. This T /T 

ratio is consistent with previous experimental thrust data which 

implied the electrothermal thrust fraction of the quasi-steady 

MPD arcjet was constant and equal to 20% over a wide range of 
125 

conditions . 

Several implications of these thrust data may now be dis- 
cussed. Since each flow region accelerates three grams per 
second argon, the mass-averaged velocities corresponding to 
the inner and outer flow thrusts are. 
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-I 

u 

= T 1 * / m I = 

12.5 km/sec 

-0 

u 

_o , . 

= T / m o = 

6.2 km/sec 


Comparison between these mass averaged velocities and the 
anode orifice exit velocities of the inner and outer flows 
suggests that approximately 28% of the arcjet's thrust is 
derived from the flow expansion in the downstream exhaust. 

In addition, the mass averaged velocity u 1 is in close 
agreement with the terminal exhaust velocity of 12.3 ^ 1.1 
km/sec measured in the plume of the inner flow. This veloc- 
ity also equals the maximum velocity obtainable from the 1-D, 
nozzle-like expansion of the inner flow discussed previously. 
From this correspondence, it may be inferred that thrust re- 
covery in the downstream exhaust regions of the discharge re- 
sults from the experimentally observed electromagnetic con- 
finement of the flow. 

Finally, the specific impulse of the two propellant 
streams may be calculated. 

I 1 = T 1 / m T g = 1280 sec. 

sp ' 1 3 

I° p = T°/ m Q g = 640 sec. 

The inner and outer flow regions delineated by the two-flow 
plasma acceleration model may be further characterized by the 
disparity between their specific impulses. The inner flow 
region generates thrust with a specific impulse which con- 
siderably exceeds the AVCO limiting I of 890 seconds for the 
quasi-steady acceleration of argon propellant. In contrast, 
the specific impulse associated with the outer propellant flow 
falls beneath this limiting, value . On combining these two 
flows, the overall specific impulse of the thruster, equal to 
the total thrust divided by g times the total mass flow, 
becomes , 
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I = (T 1 + T°) / g (ftj + * 0 ) = 960 sec. 

which, although lower than I* , still exceeds the minimum 
power model's fundamental limit of 890 seconds. 

The implications of the two— flow model and the above 
results, with respect to improving the performance of the 
quasi-steady MPD thruster, are clear. Namely, higher 
accelerator thrust densities and specific impulses should be 
achieved if the relative importance of the outer flow region 
is reduced. This could seem to be accomplished by reducing 
the anode and chamber radii such that the outer flow dis- 
charge chamber region is physically eliminated. If as a re- 
sult of the chamber geometry change, the inner flow properties 
remain the same, the specific impulse of the accelerator 
would increase from 960 to «i 1280 seconds. It must be re- 
called, however, that previous variations in the mass injec- 
tion geometry resulted in increased terminal voltages and 
erratic arcjet behavior when the outer anode regions were 
deprived of injected mass. Furthermore, power density 
and luminosity measurements imply ionization phenomena are 
important in the vicinity of the anode. These data all tend 
to suggest the function of the outer flow may be primarily 
related to current conduction processes in the anode region. 
Since these phenomena and their connection and influence upon 
the plasma acceleration process are not well understood at 
this time, it is not clear that this flow region can be elimi- 
nated without inducing other significant consequences through- 
out the remainder of the discharge. Improved MPD thruster per- 
formance realized by enhancing the inner flow characteristics 
of the discharge at the expense of the outer flow, thus remains 
to be shown. 
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B. Cathode Studies (Rudolph) 

In previous photographic and diagnostic studies of the 

i'lPD discharge, it was observed that an intensely luminous 

plasma was formed in the region immediately downstream of the 
141 

cathode tip. It was postulated that this plasma was a 

result of the high axial current densities at the cathode tip 
which caused a large, radially inward "pumping" force, in- 
creasing the number densities and temperatures in this region. 

The role of this intense plasma in the overall acceleration 
process is unclear. According to the self-field thrust equa- 
tion, a positive thrust increment is generated within this 
region by the plasma pressure acting against the conical cath- 
ode surface. On the other hand, the high temperature and 
pressure in this volume elevates the propellant to excited 
states which in turn leads to radiation and frozen flow losses. 

In order to clarify the role of this cathode tip plasma, 
a much longer cathode was installed in the MPD apparatus. It 
was anticipated that with this cathode, the current would be 
distributed such that the axial current component causing the 
radial pumping would be eliminated, and the operation of an 
accelerator of a purely "sweeping" type could be studied. 

An additional advantage of this long cathode configuration 

arises from its relation to a series of experiments recently 

155 * 

conducted by Boyle. In this work, the discharge current J 

at which the arc reached a seeming limitation on quiescent opera- 
tion (as indicated by a rapidly fluctuating or "hashy" terminal 

voltage V*) was found to be a linear function of the cathode 

2 

surface area over a range of areas from 5 to 40 cm . Throughout 
this investigation, the discharge current (typically 10 to 20 kA) 
was determined to attach over the entire cathode surface, indi- 
cating a decreasing surface current density for larger cath- 
odes. The new long cathode provides an opportunity to examine 
the "critical" current for a cathode with a surface area much 
larger than those previously tested. 
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The cathode configuration initially used in this study 
is shown in Fig. 19a, it consists of a 13.6 cm long, 1.9 cm 
diameter cylinder with a conical tip, 2.54 cm in length. The 
diameter is the same as previous cathodes in order to maintain 
the same geometric constant in the self-field thrust equa- 
tion. 61 The cathode is made of type . 304 stainless steel 
instead of the usual thoriated tungsten due to the former's 
availability and relative ease of machining. The rest of the 
discharge configuration was kept the same as Boyle's to facili- 
tate a comparison with his data. The mass flow was 6 g/sec of 
argon propellant. 

A current-voltage characteristic for this configuration 
is shown in Fig. 20. Two interesting features of this graph 
are the discharge voltage of 115 volts at a current of 16 kA, 
and the voltage, V*, at the initiation of the “hashy" voltage 
of 216 volts. These values can be immediately compared to 
Boyle's data. 

A graph of discharge terminal voltage versus cathode 
surface area is shown in Fig. 21 for all previously tested 

tungsten cathodes, which were limited to surface areas of 

2 2 

40 cm . The points at 75 cm represent the present 13.6 cm 
long stainless steel cathode. The dashed lines show that the 
present data are not inconsistent with an extrapolation of the 
Boyle data. 

In Fig. 22, the onret of the "hashy" voltage is graphed 
versus cathode surface area. The value for the 13.6-cm long 
cathode is also shown and again agrees with a reasonable extra- 
polation of previous data. Both graphs indicate that the cath- 
ode emission processes may be independent of cathode material. 


i 
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o) 13.6 cm CATHODE 



b) 2S.3 cm CATHODE 


CATHODE GEOMETRIES 


FIGURE 19 
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FIGURE 21 
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CRITICAL VOLTAGE vs CATHODE AREA 


FIGURE 22 

AP 25-5074 
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The graph of the critical current versus the cathode area 

2 

is shown in Fig. 23. Again, the data for areas below 40 cm 

2 

are from Boyle while the bar at 75 cm represents the cathode 
used in the present study. If the cathode material is un- 
important, the new point indicates than as the cathode area 
increases, the previously rising critical current approaches 
an asymptotic value of about 25 kA. Further data are required 
to verify this trend. 

In order to determine whether the previously mentioned 
electromagnetic compression region off the cathode tip was 
eliminated using this larger cathode, the local magnetic field 
distribution was measured for a total current of 16 kA. The 
resulting current contour map is shown in Fig. 24. It is 
noted that at the cathode tip, there is a sizeable fraction 
of the total current (approximately 4 kA) that has a signifi- 
cant axial component. This indicates that the "pumping" force 
has not been completely eliminated in this downstream region. 
Current conduction upstream of the 6 kA line is virtually radial, 
however, and seems to indicate that a completely radial current 
distribution is possible if the cathode attachment is allowed 
to spread out over a larger area. 

The total current attachment downstream of a point on the 
cathode surface versus the axial coordinate of that point is 
shown in Fig. 25. As can be seen, a straight line can be drawn 
through the data on the cylindrical part of the cathode, which 
is consistent with a uniform attachment. 

The 2 kA attachment on the conical cathode tip and the large 
axial component of the current downstream of the 4 kA stream- 
line both tend to indicate that a longer cathode is needed to 
provide a current attachment free of cathode tip influence and 
to eliminate the radial j x B force. Consequently, a larger 
stainless steel cathode, 26.4 cm in length, has been built and 
is shown in Fig. 19b. This cathode is presently undergoing the 
same tests as the previous 13.6 cm cathode. 
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ENCLOSED CURRENT CONTOURS 


FIGURE 24 
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C. Terminal Voltage Fluctuations in Quasi-Steady MPD > 

Discharges (Boyle) \ 

The onset of fluctuations in the terminal arc voltage of 
the quasi-steady MPD accelerator has been reported previously. 

Originally this behavior was interpreted as indicative of a 1 

fundamental limitation imposed upon the quasi-steady accelera- 
tion of a plasma flow. This arcjet model proposed nominal 

MPD discharge operation was restricted to arc currents and 

2 - 

injected mass flow rates for which the parameter J /m, the arc 
current squared divided by the propellant flow rate, was equal 
to a critical constant dependent only upon accelerator geometry 
and propellant species. Empirically this critical J /m value 
equalled (J ) /m, where J is the arc current at which the onset 
of terminal voltage fluctuations is observed. 

* 2 

Recent experiments however suggest the critical (J ) /m 
parameter is neither constant nor unique. L ' > ~ > In fact, J is 
found to be functionally dependent upon both cathode surface area 
and mass injection geometry for a given propellant flow rate. 

These data also suggest the phenomena responsible for the voltage 
fluctuations occur in the vicinity of the cathode electrode. ' 

As part of a preliminary investigation to determine the 
nature and significance of these phenomena, local potential mea- 
surements about the electrode regions of the discharge were made. 

The accelerator configuration for these measurements included a 
7.62-cm-long tungsten cathode, 50:50 flow division between cath- 
ode base annulus and outer injector holes and 6 g/sec total mass 
flow rate. Two floating potential probes are positioned close to 
the central cathode and outer anode lip. The potential drop be- 
tween the cathode and the potential probe several millimeters 
off the cathode surface is measured differentially and labelled 
V^. The floating potential several millimeters off the anode 
lip is simultaneously recorded with respect to the anode ground 
and is labelled V ft . The potential difference across the plasma 




flow region, V p , may be expressed as: 



where V is the terminal arc voltage. This experimental 
arrangement is schematically reproduced in Fig. 26. 

Figure 27 compares the oscillogram signatures of V, 

V A and V c for an arc current greater than the onset current 
J*. The cathode voltage signal V c clearly reproduces the 
onset time, frequency, and amplitude of the terminal voltage 
fluctuation. This 1:1 correspondence implies the fluctuations 
in the terminal voltage signal are localized about the 
vicinity of the cathode. 

The onset of terminal voltage fluctuations is further 
characterized by a sudden increase in the local electric field 
about the cathode. Figure 28 presents V c , and the terminal 
voltage V as a function of current. The current at which the 
cathode voltage V c sharply increases corresponds exactly with 
J* determined from the terminal voltage-current characteristic. 

The anode voltage V A also exhibits a similar sharp increase in 
magnitude although at a lower current of 8 kA. Coincidentally, 
this value approximately equals the ion current which would be 
generated by the complete single ionization of the outer 3 g/sec 
argon mass flow rate. No such sharp voltage increases are evi- 
dent in the terminal voltage and V p characteristics. 

It is interesting to note that the terminal voltages at 
which the sudden increase in cathode potential and discharge 
instability occur, V , are insensitive to the cathode surface 
area. Preliminary stainless steel large cathode results suggest 
that V is also insensitive to cathode material (section III-B) . 
Referring back to Fig. 22 , the common voltage associated with the 
onset of the voltage fluctuations is 180 volts ± 30%. Thus it 
would appear that for arc conditions such that the terminal voltage 
exceeds V*, the local electric field about the cathode sharply 
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increases and becomes oscillatory. It is tempting to specu- 
late that this critical voltage and the changes in the field 
properties about the cathode reflect a change in the local 
current emission processes. However, since the electron 
emission and current conduction processes in the cathode 
region of an arc discharge are extremely complex and perti- 
nent experimental data do not presently exist, further clari- 
fication of this suggestion must await future work. 

t * * 

Finally , there is an indication that V phenomena 
may ultimately influence the performance characteristics of 
the quasi-steady MPD arc jet. Figure 29 illustrates the 
functional relationship between the thermal efficiency W(-h 
and arc current for the above accelerator configuration. 
Assuming that JV C and JV A represent all the power that is 
deposited into the cathode and anode respectively, then the 
thermal arcjet efficiency neglecting radiation is defined as 


A 


Hi 




Physically, ^ represents the fraction of input energy 
which appears as plasma enthalpy during the arc discharge. 

From a performance point of view, ^ equals the maximum 
attainable arcjet thrust efficiency for the ideal case wherein 

all the plasma enthalpy is converted into useful thrust. For 

* 

arc currents less than J , ^ monotonically increases. Above 

J , however, the thermal efficiency assumes a constant value of 
approximately 6f%. Whether or not this trend repiesents a 
fundamental limitation imposed by electrode phenomena upon 
quasi-steady arcjet performance deserves further attention. 
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IV. PLASMADYNAMIC LASER STUDIES 

MPD Discharge in a Laser Cavity (Dutt) 

This section describes progress on a set of experiments 
intended to determine directly whether lasing can be sustained 
in the discharge and exhaust regions of an MPD accelerator. A 
direct evaluation of lasing is possible by placing the region 
of plasma under consideration within a resonant optical cavity. 
In order to minimize the variation of plasma properties along 
the optical axis of the resonant cavity, a two-dimensional 
discharge geometry was considered desirable. 

155 

The previous semi-annual report described the important 
characteristics of the cavity. Progress described in this re- 
port includes a) the alignment of the cavity mirrors to form 
a resonant cavity; b^ the measured distributions of radiance 
and current in an ambient two-dimensional discharge; and c) the 
design and calibration of a mass injection system capable of 
providing a short risetime, uniformly distributed mass pulse 
in the two-dimensional geometry. 

Ali g nm ent of the Cavity Mirrors . 

Two concave mirrors 1.5 cm in diameter with a radius of 
curvature of 67.5 cm constituted the resonant cavity. Their 
reflectivities wer 97% and 99.7% in the wavelength range 
430 nm to 530 nm, which includes many of the strong emission 
lines of ionized argon, krypton, xenon and many other elements. 

An externally mounted helium-neon laser, wavelength 
632.8 nm, was used to align initially the cavity. To ensure 
that the axes of the cavity and He-Ne laser were coincident 
and parallel to the supporting I-beam, adjustable irises 
mounted on an optical benph were used. For the initial align- 
ment, the multiple reflections of the laser beam on the inner 


i 



-63 


reflecting surfaces of the mirrors were made to coincide. 

The He-Ne laser was then removed. Subsequently, a commercial 
argon ion laser, with its mirrors removed, was placed within 
the cavity (Fig. 30a) . The laser tube was lined up with the 
cavity axis and the discharge initiated. Laser output was 
observed from the mirror of lower reflectivity. Minor adjust- 
ments in the mirror positions were made to obtain the brightest 
output. Each mirror could be rotated approximately 0.25° about 
a diameter before lasing disappeared completely. The mirrors 
were easily adjustable to within 0.02°. The argon ion laser 
tube was then removed. 

Following alignment, the cavity can now be used to in- 
vestigate whether conditions suitable for lasing (described 

155 

in the previous semi-annual report ) exist in any medium 
within the cavity. 

Two-dimensional MPD Discharge 

Experimental arrangement . The coaxial geometry predomi- 
nant in MPD arc work was considered inadequate to provide a 
lasing medium because of the large gradients of plasma prop- 
erties along any line of sight perpendicular to the MPD axis. 
Specifically, the gain of an amplifying region could be masked 
by other absorbing regions elsewhere along the line of sight. 

In addition, a longer active medium is prefeiuble because the 
overall gain (or loss) within the medium is correspondingly 
increased. This would also provide a more accurate determina- 
tion of the gain per unit length. 

To realize these advantages, a two-dimensional version cf 
the MPD arc was designed and constructed with electrodes con- 
siderably wider than any previously used in this laboratory. 
This was intended to provide a long uniform plasma region with 
variations in the plasma properties confined to relatively 
small end regions. 
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o) RESONANT CAVITY 



b) DISCHARGE APPARATUS 
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The proposed discharge geometry was briefly described 

155 

in the previous semi-annual report. It consists of three 

aluminum parallel plates, 45-cm wide by 0.9b-cm thick as 
shown in Fig. 30b. The upper and lower plates are anodes 
which extend 10 cm out of the backplate, while the center 
cathode plate extends only 5 cm from the backplate. The 
plates are tapered at their rear ends to fit into a 20-cm 
diameter Plexiglas cylinder which passes through the circular 
back access port in the large Plexiglas tank. Electrical con- 
tact to the plates from outside the tank is made by brass studs 
passing through the Plexiglas cylinder. Copper plates connect 
the brass studs to the capacitive pulse-forming network via 
a gas triggered switch. 

When the assembly is in place within the vacuum tank, 
the electrode region lies in ..he line of sight of an opposing 
pair of tank windows. The resonant cavity may thus be placed 
with the mirrors outside the windows to contain various por- 
tions of the discharge region. 

In order to provide a suitable lasing medium, the dis- 
charge must satisfy several criteria. One important require- 
ment is that the components of the current to the two anodes 
be equal for the discharge to be horizontally symmetric. It 
should be noted that the two anodes are not physically con- 
nected inside the vacuum tank. Second, the current distribu- 
tion must be diffuse and approximately uniform across the 
entire width of the electrodes with relatively small fringing 
effects. Finally, not only must the current be uniformly 
distributed but mass must also be injected uniformly across 
the electrode width in order to guarantee uniform plasma 
properties in the cavity. 

In order to determine whether these criteria are satis- 
fied for our discharge configuration, several experiments 
were conducted in an ambient gas discharge. This technique 
has been found useful in the past for examining the feasi- 
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bility of a particular experimental configuration without 
the added complexity of a mass injection system. 

Ambient Gas Discheirqe Operation . For these initial tests, 
the vacuum tank was filled with argon to pressures varying be- 
tween 25 fX and 300 /a of mercury. This is also the range of fill 
pressures for conventional argon ion lasers. During a discharge, 
the terminal voltage was measured using a Tektronix voltage 
probe. The current distribution between the two anodes was mea- 
sured by passive integration of separate Rogowski coil signals. 

Oscillogram traces of current and voltage are shown in 
Fig. 31 for a total current of 65 kA supplied for 200 fxsec. 

The two uppermost traces represent the currents and to 
the upper and lower anodes, respectively. The other traces 
are voltage records for various ambient pressures. 

Current measurement over a variety of operating conditions 
show the current division to be quasi-steady, and in all cases 
and J 2 differ by less than 5%. The device thus satisfies 
the important condition of equal division of current between 
the two anodes. 

The voltage traces exhibit the familiar drop to a minimum 
after discharge initiation, followed by a rise to a quasi- 
steady voltage level. These features are usually associated 
with insulator ablation, in this case caused by operation in 
the ambient mode. 

Photographic Study . Photographs of the discharge were 
taken for a limited range of conditions as indicated in the 
table below. 


Table II. Photographs of Ambient Discharge 


Figure 

Number 

Current 

Amplitude Duration 

Camera Lens 
Aperture 

Neutral Density 
F: Iter 

32b 

8.0 kA 

0.92 msec 

f/32 

50% 

32c 

15.6 

0.92 

f/32 

50 

32d 

15.6 

0.92 

f/H 

10 

32e 

65.0 

0.20 

f/16 

10 

33b 

15.6 

0.92 

f/11 

10 
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At 8 kA the luminosity is confined to one— half of the 
discharge width (Fig. 32b) indicating that this current is 
inadequate to drive the 45-cm wide discharge. At a higher cur- 
rent of 15.6 kA the discharge appears to be distributed more 
evenly but it is still much brighter on one side. Fig. 32c. A 
different photograph of this same discharge in Fig. 32d taken at 
a different combination of lens aperture and neutral density 
filter shows it to be more uniform. When the current was in- 
creased substantially to 65 kA, the discharge appeared more 
laterally symmetric but exhibited bright areas near the ends 
of the electrodes. Fig. 32e. 

A side view at 15.6 kA (Fig. 33b) shows the discharge to 
be horizontally symmetric. The bright radiance close to the 
ends of the electrodes originates on the surfaces of the in- 
sulating sidewalls, as can be seen by comparing this photograph 
with the no-discharge photograph shown in Fig. 33a. 

The more uniform distribution of luminosity which is ob- 
served at the higher currents may be taken as an indication of 
a relatively uniform distribution of current and plasma para- 
meters, but it is not in itself sufficient evidence of uniform- 
ity. It does suggest however, that the plasma properties are 
more uniform at higher current levels. This favors the choice 
of high current levels for a finer and more extensive mapping 
of the current distribution in the discharge. 

Current Distribution . Measurements of magnetic field may 
be used to determine not only the uniformity of current distri- 
bution across the discharge but also to evaluate the current 
densities within the discharge regions. 

To facilitate the reduction of magnetic field data a 
coordinate system shown in Fig. 30b was used with the origin 
centrally located on the cathode tip. Three similar probes 
were used, instead of one, to speed up the probing of such an 
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extensive discharge region. The probes were mounted on a 
movable platform within the vacuum tank with their stems orien- 
ted parallel to each other and pointing in the negative z 
direction. The axes of the probe coils were in the x direc- 
tion so that the component B x of the magnetic field would be 
measured. This was expected to be the dominant component for 
the given orientation of electrodes. Probes #1 and 2 were on 
the same horizontal plane 21.6 cm apart while probe #3 was 
5.7 cm directly below probe #2 (Fig. 34a) . The probe coils 
were all at the same z value. 

The magnetic field B x was measured at a few locations in 
the discharge at the intermediate current level, 15.6 kA. Os- 
cillograms reproduced in the left column of Fig. 34b represent 
magnetic probe signals at locations x = -21.6, 0 and 4-21.5 cm 
on the line y = 4 cm, Z=-3.8 cm. Corresponding oscillograms for 
the 65 kA case are shown in the right hand column of Fig. 34b. 
The traces indicate that at the lower current, the field com- 
ponent B x and consequently the current distribution remain un- 
steady for a large part of the 0.9 msec pulse. At the higher 
current, 65 kA, the discharge becomes quasi-steady in about 
120 psec. 

The distribution of the quasi-steady values of B x in the 
x direction is shown in Fig. 35 for probing along three lines 
at the y=2 cm plane. These locations correspond to the insulat 
ing backplate Z= 3.8 cm, the cathode front edge, Z = 0, and the 
anode front edge, Z= 5.1 cm, respectively. 

The most important conclusion from these data is that the 
discharge is uniform across the entire 45-cm width of the 
electrodes at locations close to the insulating backplate. Fur- 
ther away from the backplate the profile is less uniform with a 
B x of one-half its maximum value recorded at the edges of the 
electrodes. The reason for this decrease in the field may be 
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that the insulating sidewalls, which are partially cut away 
(see Fig. 30b), permit the current to fringe near the ends of 
the electrodes. The lateral fringing of current contours is 
substantiated by measurements of B which have values comparable 
to B x at the edges of the discharge, downstream from the back- 
plate. It is expected that this small departure from uniform- 
ity will be corrected by using more extensive sidewalls. 

To evaluate the vector current density j, all the com- 
ponents of the magnetic field within the discharge must be 
known. However, in these initial experiments only B x , the 
dominant component of the magnetic field, was measured in some 
detail over the entire region of the discharge. For a total 
current of 65 kA these measurements lead to an average value of 
j = 93 A/cm between the cathode and anode in the plane contain- 
ing the anode front edges about the symmetry, or 2 , axis. It 
the other components of B had been included, j would be some- 
what larger. It should be noted that the value of j is in 
excess of the typical threshold current density, 50 A/ cm , for 
conventional argon ion laser operation, although the signifi- 
cance of this comparison is not yet known for this discharge con- 
figuration. 

The experiments with the ambient argon fill discharge opera- 
tion show that a fairly uniform discharge of sufficiently high 
current density can be sustained. However, the voltage measure- 
ments indicate the possible presence of substantial ablation 
products in the discharge. In order to operate the discharge in 
a pure argon plasma, a means of injecting uniformly a sufficient 
mass flow of argon gas into the discharge must be provided. The 
design of such a system is presented in the following section. 

Mass Injection System 

Design, The mass injection' system should provide a uniform- 
ly distributed quantity of gas into the discharge region at a 
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rate that enables a cniasi-steady discharge to be maintained 
with minimal ablation. From previous experiments with a 15-cm 
wide parallel plate accelerator, the mass flow required for the 
highest current anticipated (^80 kA) may be greater than 
100 g/sec. 128 

The risetime of the mass pulse is also crucial in that the 
mass injected prior to discharge initiation raises the back- 
ground pressure in the vacuum tank. To the extent that this 
background pressure can be compared to the initial fill pressure 
of conventional argon ion lasers, the riset.ime takes on added 
significance. Conventional argon ion lasers operate with a cold 
gas filling pressure in the range from 25 to 300 p. In com- 
parison, a 100 g/sec mass pulse with a 10 msec risetime will 
raise the average pressure in the present vacuum tank to approxi- 
mately 400 p. Although a strict comparison between the two 
pressures is not warranted because of the differences between the 
conventional laser and the present dynamic system, it is never- 
theless apparent that the risetime should be kept as small as 
possible . 

The mass injection system is shown in Fig. 36a. A standard 
gas cylinder was positioned immediately outside the vacuum tank 
and connected with a tee arrangement to the two solenoid valves 
fixed to the side of the ’ nsulating backplate. The valves are 
Skinner type V52 with a nominal diameter of 0.8 cm. Holes 
drilled through the side of the insulating backplate connect the 
valves to two plena which run across the backplate. Each plenum 
is a 0.95-cm square channel machined out of the Plexiglas back- 
plate. The plena are closed by two rectangular plates (0.63-cm 
thick) which cover the i nterolectrode region. There are nine 
equidistant holes of 0.159-cm diameter on each plate directly 
in front of the plenum. 

The solenoid valves are opened by discharging an 850 pF 
capacitor through them. Earlier experiments have shown that the 
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two valves open almost simultaneously, delivering gas into the 
two plena. The nine holes in each plenum allow the gas to 
escape and fill the discharge chamber. 

2 

The effective choking area of each valve is 0.196 cm . The 

2 

total area of the 9 holes is 0.181 cm . Thus, irrespective of 
the discharge coefficient of the holes, the flow when steady 

will choke at the nine exit holes. The advantage of such an 
arrangement is that when the discharge commences the flow adjust- 
ment time is small, thus allowing quasi-steady operation to be 
reached in a shorter time. 

Performance . A fast ionization gauge (Varian Millitorr) 
was placed at various locations within the tank in order to 
determine the uniformity of the cold flow. The time for the 
flow to reach a quasi-steady state was measured in addition to 
the distribution of pressure once the steady state was reached. 

A representative oscillogram of the ionization gauge out- 
put is shown in Fig. 36b. The pressure at the anode region 
reaches a steady V'lue in about 2 msec. A profile of steady 
state pressure distribution at y = 3 cm, Z = 5.1 cm is shown 
in Fic . 37. The peak- to- trough variation is j+157 about the 
mean, while the decrease in the average pressure from the left 
to right is 8% consistent with the valve location on the 
left or negative x end. These data show that the flow reaches 
a steady state quickly and is reasonably uniform in the x 
di rcction. 

Calibration of the Mass Flow Rate. The mass flow was cal- 

119146 

ibrated as described in earlier semi-annual reports. ' 

A Hewlett Packard Pulse Generator was used to generate a square 
wave lasting 2.5 to 10.1 msec. The square wave was differen- 
tiated, with the resultant positive and negative spikes used 
to tailor the voltage pulse to the valves. The positive spike 
allows the previously charged 850 /iF capacitor to discharge 
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in to the valves, opening them after a small delay. The nega- 
tive spike shorts out the capacitor allowing the valves to 
close. With the vacuum pumps isolated, the pressure rise in 
the tank due to the injected mass pulse is measured using a 
McLeod gauge. 

By altering the length of the square wave and thus the 

Lime interval between the spikes, the mass of gas admitted into 

the vacuum chamber can be varied. A minimum pulse length t . 

exists below which the valve never fully opens and thus the 

flow rate through the valve does not reach the choked flow rate 

associated with the particular stagnation pressure. For pulse 

lengths greater than the minimum, a fixed mass flow is obtained, 

and the tank pressure increases linearly with t — t . . The 

mm 

slope of the tank pressure rise vs. pulse length for time 
greater than the minimum thus yields the mass flow rate. 

A typical calibration curve gives an average pressure rise 

5 2 

of 1.81 /t/ msec for a stagnation pressure of 1.34(10 ) N/m (5 psig) 

3 

Using the tank volume of 1.22 m and an assumed ambient temper- 
ature of 295 °K, the mass flow is calculated to be 4.74 g/sec. 
Similar measurements at double the absolute pressure give twice 
the mass flow verifying the choked flow assumption. 

Summary . 

The results and conclusions of the experiments described 
in this report may be summarized as follows. 

Alignment of the constructed optical cavity with a com- 
mercial argon ion laser source is readily accomplished. This 
cavity is now being adapted for first experiments with the 
large Plexiglas vacuum tank. 

The two-dimensional discharge geometry was constructed 
and tested in ambient argon- at pressures around 200 / a . Rogowsky 
coil measurements showed that in all cases the current was dis- 
tributed evenly between the two anodes. In addition, for 
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suf ficiently large currents the discharge distributed itself 
to provide a uniform luminosity and magnetic field along the 
width 01 the electrodes. The magnetic field distribution was 
less uniform away from the insulating backplate, probably 
because restraining sidewalls were absent there. Large side- 
walls are presently under construction. 

For a total current of 65 kA the current density in the 
middle of the discharge was estimated to be twice the threshold 
current density required for argon ion laser operation in thin 
tubes. About half of the total curren was found to attach on 
the outer surfaces of the anodes. This suggests that the cur- 
rent density could be increased in some parts of the discharge 
by insulating the outer surfaces of the anodes. 

As was expected, ablation of the Plexiglas insulators oc- 
curred during the discharge, especially at the larger currents. 
These insulators will be replaced or covered with some refrac- 
tory material such as quartz or boron nitride. This will re- 
duce the contamination and allow a clearer interpretation of 
the observed phenomena. 

A mass injection system was designed and constructed. 
Measurement of cold flow pressure distribution across the dis- 
charge region showed it to be reasonably uniform. The mass 
outflow was experimentally determined. Sufficiently high mass 
flows can be obtained at stagnation pressures of several at- 
mospheres with the possible addition of two other valves. 

These results show that a quasi-steady uniform two-dimen- 
sior.al discharge can be obtained over wide plates, accompanied 
by uniform gas injection at desired levels. Experiments can 
thus be performed to determine whether lasing is possible in 
the medium under various conditions. 


i 
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V. HOLLOW CATHODE STUDIES (Krishnan) 

Advantage was taken of the recent move of the Electric 
Propulsion Laboratory to redesign the pulse-forming network 
to generate longer current pulses. The new pulse lengths 
are 500 p.s ec in the parallel and 2 msec in the series con- 
figurations as compared with 200 ^sec and 800 jisec, respec- 
tively, from the previous arrangement. The discharge cham- 
ber and its supporting equipment were also rearranged to 
allow greater accessibility for photographic and spectro- 
scopic study of the hollow cathode discharge. 

In earlier studies the distributions of current and 
potential within and about hollow cathodes of different con- 
figurations, designated as HC-1 to HC-VIII, were determined 

in order to establish some of the characteristics of typical 

151 155 

hollow cathode operation. ' These studies have been 

continued with two additional hollow cathodes, designated 
HC-X and HG-XI , in order to determine the distributions of 
current and potential which distinguish typical hollow cathode 
operation in a high current discharge. For reasons of ex- 
pediency, the most recent two hollow cathodes HC-X and HC-XI 
were manufactured of stainless steel instead of tungsten as 
were all the previous configurations. Therefore, in order 
to compare the experimental results with those acquired 
previously, one additional stainless steel cathode, HC-IX, 
of the: same dimensions as the tungsten HC-VIII, was manu- 
factured to record any differences in its electrical char- 
acteristics, as well as to check out the operation of the 
improved facility in the relocated laboratory. 
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This report describes in its first section the modifica- 
tions and improvements of the hollow cathode facility. The 
second section presents the results of studies with the two 
additional new hollow cathode configurations. These studies 
bring to a conclusion the first phase of the investigation 
of the characteristics of hollow cathode operation in a high 
current discharge. 

Expand ed Hollow Cathode Discharge Facility 

The major change in the hollow cathode facility is an 
extension of the pulse-forming network to provide longer cur- 
rent pulses at the same current levels as before. The previ- 
ous pulse-forming network provided approximately 4 to 20 kA 
for 200 psec when arranged in the parallel configuration or 

1 to 5 kA for 800 psec when in the series configuration. The 

143 

old facility, which was described m a previous report, 
consisted of 4 rows of individual transmission lines of 9 
capacitors each, giving the capacitor bank a total capacitance 
of 940 pF. 

The new facility consists of 4 individual transmission 
lines, each assembled from 21 equal sections of 6.6 pH series 
inductors and 27.1 pF shunt capacitors as shown in Fig. 38. 

The total capacitance of this new bank is 2.35 mF. With this 
larger network the line segments can be arranged in parallel to 
provide approximately 4 to 20 kA for 500 psec or in series to 
provide 1 to 5 kA for 2 msec, thus yielding current pulses 2-1/2 
times as long at the same current levels as the older and smaller 
capacitor bank. A photograph of the capacitor bank is shown in 
Fig. 39. 

Voltage reversal across any of the capacitors in the pulse 
forming network must be avoided to prevent damage. In addition, 
even a low level of reversed current is undersirable since it 
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changes the polarity of the electrodes of the MPD apparatus. 

'The resulting change in the distribution of radiance from the 
discharge completely masks the actual quasi-steady state of the 
discharge in photography and spectroscopy unless high speed 
resolution techniques are used. 

To prevent a voltage reversal in the network due to the 
mismatch between the network impedance Z N and the 10 mfl ac- 
celerator impedance Z^, a resistor R must be connected in 
series between the network and the accelerator as shown in 

Fig. 38 . If the value of the resistor R = Z. - Z,. , then the load 

N A 

formed by the accelerator and the resistor is effectively matched 
to the network and no appreciable voltage reversal occurs when 
the network is discharged across the load. 

The series and parallel configurations of the ladder net- 
work use '■* ectrolytic load resistors of 410 mo and 122 mn re- 
spectively. These resistors, shown in the lover left of Fig. 

39, consist of 2 copper plates 0.08 cm thick, 33 cm x 63 cm, 
with a separation of 5 cm suspended in a polyethylene tub filled 
with 56 liters of cupric sulphate electrolyte. Besides their 
simplicity and low initial cost, these resistors are particular- 
ly advantageous for this applicat; on due to the very large heat 
capacity of the liquid. It is estimated that the temperature of 
the solution rises by less than 0.50 °C during the 20 kA x 500 
psec pulse. 

The discharge chamber, glass vacuum tank and associated 
electronics are positioned on a table in the middle of a separate 
room as shown in Fig. 40, thus giving easy access from every 
direction for spectroscopy and photography. Furthermore, the 
connection between the gas triggered switch of the discharge ap~ 
paracus and the wall outlet of the power supply network is made 
by ten flexible low inductance RG14-A/U coaxial cables, shown 
in Fig. 40, of sufficient lengthto allow the entire accelerator 
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assembly, including the table, to be rotated by up to 93° 
when needed for diagnostic study of the discharge. 

Two typical oscillograms of the current pulse from the 
new pulse forming network are displayed in Fig. 41. The upper 
trace, Fig. 41a, shows the 6 kA x 500 psec pulse when the 
parallel netwr rk ir initially charged to 1.7 kV. The droop 
cf the 2.9 kA x 2 msec current pulse from the series network 
in Fig. 41b is due to the relatively short time constant, 

= 4 usee, of the R-C integrator which was used for recording 
the signal from the Rogowski coil. 

E xperimental 

Previous semi-annual reports contained the results of 
detailed experimental investigations of eight difterent in- 
sulated hollow cathode configurations designated HC-I to 
HC- VIII . Before describing recent work with tne new hollow 
r ’ thode configurations, designated HC-IX, X and XI and sh iwn 
in Fig. 42, it may be worthwhile to recall the important re- 
sults of the studies of cathodes I to VIII in order to under- 
stand the need for further alterations of the cathode geometry. 

These features are: a) Electron . number densities in excess of 
17 - 3 

10 cm prevail in the cavity of the hollow cathode; b) A field 

free region with a weak potential gradient of less than 10 v/cm 

exists in the same region of the cavity where current densities 

2 

of typically 1500 A/cm are extracted from the wall; c) The 
observed correlation of axial electric field E in the in- 
sulated channel with area A of the insulated channel 

E (V/cm) - 6.15/A (cm 2 ) (V-l) 

lends credence to the assumption of constant electrical conduc- 
tivity in the channel of 10 mho/m which corresponds to an elec- 
tron temperature of 3 to 4 electron volts at any substantial 
ionization level. 


y. 


- 68 - 


— — j |- — 100/u.sec 


1C • 1147 



— H |* — 500 sec 


IC-126S 



b) SERIES 


CURRENT SIGNATURES 


6kA/DIV 


2.4 kA/DIV 


FIGURE 41 

AP25P-542 











-90- 


While these are interesting observations they do not yet 
establish a complete set of criteria for hollow cathode opera- 
tion. In an effort to identify the principal characteristics 
of hollow cathode discharges, a few additional hollow electrode 
configurations have been constructed and their performance 
examined in the laboratory. 

Effect of Cathode Material . Figures 42a and 42b show the 
previous tungsten hollow cathode HC-V1II and the new stainless 
steel cathode HC-IX of identical cavity dimensions. For this 
latter electrode, the profiles of enclosed current characteris- 
tic, were determined. As in previous studies, a current of 7 kA 
and argon mass flow of 4 g/sec entirely through the cathode 
cavity were maintained, but now with the extended capacitor bank 
the duration of the discharge was 0.5 msec instead of 0.2 msec. 
All experimental studies in this report were undertaken at these 
discharge conditions. 

For the survey of current distribution, the magnetic probe 
was translated at a radius of 0.8 cm. The measured enclosed 
current profiles for the tungsten electrode HC-VIII and for the 
stairless steel electrode HC-IX are compared in Fig. 43. The 
current was normalized to the 3 kA value measured at the down- 
stream cavity end at r=0.8 cm (0.15 cm from the cathode surface) 
The relative currant distribution inside both cathode cavities 
is seen to be identical, with more than 90% of the measured cur- 
rent attaching over the last 1 cm of cavity. The data for the 
steel cathode extend 2 cm downstream, exhibiting a pinch at 
0.3 cm followed by a rapid expansion of the current djwnstream. 

Axial potential profiles were determined experimentally 
with a floating Langmuir probe. The floating potential relative 
to the anode was subtracted from the terminal voltage signal 
with an accurately balanced difference amplifier, leaving the 
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floating potential relative to the cathode. The measured 
axial potential profiles for the tungsten HC-VIII and the 
steel HC-IX are compared in Fig. 43. The cavity is essential- 
ly field-free up to 1 cm from its end, the region of primary 
current attachment. Over this last 1 cm of the cavity the 
tungsten cathode shows a small average field of a few volts/ 
cm, while the steel cathode shows a slightly greater average 
field of 10 V/cm. The difference in potential increases to 
approximately 15 volts at an axial location 2 cm downstream 
of the cathode. 

These observations can be summarized in the statement that 
the change of cathode material from tungsten to steel produces 
only marginal changes in the measured enclosed current and 
potential distributions. 

In addition to the surveys of current distributions and 
floating potentials for HC-VJII and HC-IX, their voltage- 
current characteristics have also been recorded over a discharge 
current range from 2 to 16 kA. The characteristics are ex- 
hibited in Fig. 44 and are seen to coincide within experimental 
error. At approximately 10 kA both characteristics show a 
relatively abrupt change in their slopes. The nature of this 
transition is not being pursued further in these investigations 
of hollow cathodes at this time; the interested reader is re- 
ferred to an earlier discussion of such transitions in the slope 
of the characteristics, which have been observed previously for 
MPD discharges with solid cathodes . - 1 - 55 

New Hollow Cathode Co n figurations . In addition to the 
previously reported results of the studies with the hollow 
cathode configurations KC-I through HC-VIII, two new configura- 
tions designated HC-X and HC-XI, shown in Fig. 42, were manu- 
factured and their operation investigated in the laboratory. 

The cavity configurations were changed significantly in these 
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two models in order to observe differences in the characteris- 
tic features of current and potential distributions v:hich 
distinguish typical hollow cathode operation from discharges 
with solid cathodes. 

In both of these cathodes the diameter of the orifice in 
the face of the cathode and the diameter of the upstream entrance 
to the cavity, as well as the length of the electrode have been 
retained as they were in the preceding configurations, KC-VII 
to HC-IX. However, the previously cylindrical cavity now takes 
the shape of a right circular cone with a half angle of 10° in 
HC-X, while the cavity of HC-XI is formed by a shallow conical 
chamfer of 45° half-angle, fed by the upstream orifice through 
a relatively narrow, slightly expanding (2 C ) channel. hese two 
configurations may be considered as steps in the transition 
from a hollow cathode with a cylindrical cavity to a compact 
electrode with a flat face. Furthermore, the current distri- 
bution in conical cavities of increasing half-angles can be 
expected to contain decreasing components of radial current 
density if the attachment at the electrode cavity wall is 
assumed to remain normal to it. Excluding the thin cathode fall 
region, this assumption of normal attachment may be justified by 
the high conductivity of steel and tungsten, 1.3 x 10 mhos/m and 

7 

2 x 10 mhos/m respectively, which exceeds by at least two orders 
of magnitude the estimated plasma conductivity in the cavity of 

4 

10 mhos/m. Therefore, if the plasma conductivity is scalar, 
then the current attachment can indeed be assumed normal. 

With smaller radial current components a decreased axial 
component of the Lorentz body force is then available to retard 
the plasma flow in the cavity. The current and potential pro- 
files for hollow cathodes HC-X and HC-XI are therefore expected 
to exhibit a redistribution of currents and potentials as a con- 
sequence of the decreased retardation of the plasma flow. The 
magnetic field was surveyed with a small magnetic probe of 
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« 0.3 cm diameter containing 56 turns of #38 wire, and a 
single Langmuir probe was used to record the floating poten- 
tial. The results of these experiments are presented and dis- 

£ 

cussed in the following paragraphs. '? 

3 

The magnetic probe survey of cathode KC-X was carried 
out in steps of 0.1 cm, and in regions of large gradients in f 

0.05 cm steps, over an axial extent from -0.5 cm to + 1.0 cm ; 

and at radial positions from 0.5 cm to 0.95 cm. The floating | 

potential was measured on the cavity centerline. Figure 45 | 

compares enclosed current and floating potential profiles of f 

HC-X with those of HC-IX, reproduced from Fig. 43. The out- i 

line of the two electrodes is sketched in the top right of 
Fig. 45 for reference. < 

i 

Looking first at the current distribution, the origin of | 

the abscissa is the downstream end of the cavity, xvhere the 
normalization current of 3 kA was measured. Because of the { 

0. 3-cm probe diameter, the closest possible approach of the 
probe to the cavity surface was 0.15 cm. The distributions 
in Fig. 45 were obtained at this local maximum radius. It is 

seen that the distribution of current attachment in these two ■ 

electrodes is practically identical within experimental error, 1 

1. e. the change of the cylindrical cavity configuration to 

that of a 10° cone has no recognizable effect on the current ;• 

distribution. { 

■j 

The enclosed current and floating potential profiles ■] 

for HC-XI with the 45° chamfered cavity were similarly de- j 

termined and are shown in Fig. 46, again with the corresponding 
profiles of HC-IX. It should be noted in this figure that the i 

enclosed current is plotted with the coordinate Y as abscissa. 

This coordinate Y, which is indicated in the succeeding Fig. | 

47, is again the line of closest approach of the magnetic probe 

to the cavity surface (0.15 cm). The current profile of HC-IX is \ 
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normalized with 3.0 kA as before while that of HC-XI is nor- 
malized with 45 kA, which is the value measured at approxi- 
mately 0.05 cm from the origin of the coordinate Y. 

The enclosed current profile of HC-XI in Pig. 46 is seen 
to decrease approximately linearly at a slightly more rapid j 

rate than HC-IX and then to start levelling off. beyond Y = 0.4 
cm. _ .1 

The linear portion indicates equal enclosed current in- i 

crements for equal increments of the surface coordinate Y. j 

However, because the conical surface area increments are de- j 

creasing as Y increases, the local current density must be j 

increasing. The current density thus reaches a maximum and ! 

then drops as the enclosed current profile levels off at low j 

v J 

current values. A plot of the entire enclosed current pattern 

\ 

for HC-XI, Fig. 47, shows this maximum current density as a j 

bunching of the enclosed current contours on the 45° probing l 

* 

surface. The enclosed current contour pattern for HC-X, also 

presented for comparison in Pig. 47, does not show a pronounced i 

o ^ 

maximum because of its small 10 conical angle. The cause and 

significance of thi3 current density maximum on the face of the j 

conical cathode HC-X are not yet understood. ! 

•! 

Referring back to Figs. 45 and 46, the axial floating poten- t 

tial profiles relative to the cathode are presented for HC-X and j 

KC-XT together v/ith the profile for HC-IX. HC-X shows in its j 

long cavity a potential plateau of 30 volts, whj le KC-XI shows j 

a short potential plateau of a mere 4 volts in the upstream \ 

portion of the narrow cavity. More detailed potential measure- \ 

ments are required in order to establish the significance of I 

! 

these features. j 

i 

Following these tests with other geometries, one cathode si 

has been selected for a more detailed series of diagnostic tests. 1 
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Because geometry does not appear to be a major factor in ■ 

altering the key features of the current and potential dis- ] 

tributions, the cathode design selected was similar to 

HC-VIII, i.e., one with a large constant diameter cavity '• 

making it most accessible for probing. Detailed current ; I 

density and potential surveys arc now underway over a wide i i' 

range of operating conditions with a refractory insulator I 

coated version of this cathode. . i j 
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APPENDIX A: Semi-annual Statement of Expenditures 

PULSED ELECTROMAGNETIC GAS ACCELERATION 
NASA NGL 31-001-005 
1 January 1974 to 30 June 1974 


$18,337 

11,814 

2,770 

2,364 

$35,285 


II. Employee Benefits ( 21 %) 6,828 

III. Equipment 0 

IV. Materials and Services 7,678 

V. Travel 0 

VI. Tuition 2,320 

Total Direct Costs $52,111 

INDIRECT COSTS 

VII. Overhead (80%) 28,228 

Total $80,339 


DIRECT COSTS 
I. Salaries 

A. Professional 

B. Technicians 

C. Students 

D. Supporting Staff 
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